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Abstract 

A Monte Carlo simulation is performed in a cubic lattice of interacting identical Stoner-Wohlfarth 
nanoparticles. The model system is a randomly-anisotropic Heisenberg spin system with a small 
anisotropy-to-exchange ratio D/J = 3.5. The dc susceptibility, Xdc{T)i shows a Curie- Weiss-like 
transition at a temperature Tq/ J ~ 1.5, followed by a low-temperature glassy behavior manifested 
by cusps in both the zero-field-cooled and the field-cooled curves. The ac susceptibility, Xs.c{T.,oj)., 
at various frequencies, a;, shows that with decreasing temperature, a non-Arrhenius dispersive 
peak occurs at rb(c<-'), succeeded by another dispersionless peak at Tg/J Rr: 1.20 in the in-phase 
part, x'{T,uj)^ of Xa.c{T.,uj) while the out-of-phase part, ^"{T^uj), shows only one peak. A dynamic 
scaling analysis shows that the system exhibits a critical slowing-down at Tg with a quite small 
exponent zv ~ 1.65. However, no universal collapse is seen for the fully-scaled data of }("{T,uj). 
These observed behaviors are interpreted under the droplet-like hypothesis that the formation and 
development of exchange-induced correlated clusters drive ensembles of nanoparticles undergoing a 
transition from a paramagnetic order to a short-range order (SRO) at Tq, followed by a transition 
at Tg to the magnetic state in which a magnetic glassy order and a magnetic quasi-long-range order 
(QLRO) coexist. In addition, our simulation shows that the onset of the latter transition, which 
is peculiarly manifested by the dispersionless peak, occurs only for those ensembles possessing the 
anisotropy strength in the region 1.0 < D / J < 5.0. When this range is exceeded, this onset totally 
suppressed. The reason that the QLRO is prohibited in ensembles of strong random anisotropy 
may account for this phenomenon. 
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In 1973, Harris, Plischke, and Zuckermann (HPZ) [l] proposed the so-called random 
anisotropy model (RAM) for magnetism in rare-earth amorphous compounds because of the 
fact that both the Curie temperature Tq and the spontaneous magnetization Mg of the TbFe2 
amorphous alloy are substantially lower than those of its crystalline counterpart observed by 
Rhyne, Pickart, and Alperin [2i]. Since then, the model has been considered as the prototype 
for the physics of magnetic phenomena, not only in magnetic amorphous alloys but also in 
other disordered magnetic materials possessing random anisotropy s], Q, [g], among which 
the RAM has been subjected to various developments in modeling ensembles of magnetic 
nanoparticles (NPs) over the past decades tI. IsI. IqI. [io| . 

In this paper, we will apply the RAM to explore, using a Monte Carlo (MC) simulation, 
the nature of low-temperature magnetic collective behavior driven by the competition be- 
tween the exchange interaction and the random anisotropy in an NP array. The studied 
system is a simple cubic lattice of NPs, which is the simplest model among 3D ordered 



nanostructures 
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ii2i, m. Our simulation is carried out under the RAM in the sense 
that the magnetic moments of same-sized Stoner-Wohlfarth (SW) NPs [14] play a role as 
atomic spins in the original RAM model |l|. The magnetic moment can be represented as a 
vector fl = MsVfi, where Mg is the spontaneous moment and is assumed to be temperature- 
independent, V is the particle volume, and fi is the unit vector along the moment direction. 



The random anisotropy is that o 
K is the anisotropy constant 
can be expressed as 
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a single particle with a positive strength D = KV, where 
Therefore, the RAM's Hamiltonian in our study 



12 



1^ 



(1) 



which is identical to that for the original model. Here, h is the unit vector along the 
direction of the easy-axis, z is the directional vector for the Oz axis of the Cartesian system 
{Ozyz) along which the magnetic field with strength H is applied, and J is the strength of 
the nearest-neighbor exchange interaction, which is assumed to be positive (ferromagnetic). 
Notice that the terms in Eq. (1) are expressed in the convention that J, and H are 
"measured" in the dimension of energy while the scalar products inside the summation 
signs yield dimensionless values. The temperature of the system, T, is in the dimension of 
energy, too. 

The motivation for our work is associated with the recent MC work by Billoni, Cannas, 



and Tamarit (BCT's work) [15|] on the out-of-equilibrium dynamics of the model in Eq. 
([1]). Interestingly, they found such universality classes for nonequilibrium spin-glass (SG) 
dynamics 16| as the subaging scaling law and the fluctuation-dissipation relation for a 
cubic lattice of Heisenberg spins. Their work is particularly significant for RAM magnetism, 
which is still a subject of heavy controversy related to the nature of the ordered phases 
at low temperatures and its dependence on the degree of anisotropy (in other words, its 
dependence on D / J) [l7, 18]. The dependence on D/J may be summarized briefly as 



follows: in strong D/J region, the RAM seems to present a low-temperature spin-glass- 
like phase (or a speromagnetic (SPEM) phase in Coey's terminology \4\). In infinitesimally 
small anisotropy region, the order tends to be of locally noncollinear ferromagnetic order 
(also asperomagnet (ASPEM)) 4, ll9|. Weak anisotropy seems to drive the ASPEM order 



into the so-called correlated SPEM or correlated spin glass (CSG) [19|. Recently, Itakura 
{20! performed a MC simulation and depicted a phase diagram for the anisotropy dependence 
in terms of the spin-spin correlation G{r). According to that work, at low temperature, a 
quasi-long-range order should persist in coexistence with the CSG, at least up to an upper 
bound of D/J ~ 5.0 with a frozen power-law G{r) oc r^^~^, above which the spins attain 
a short-range order (SRO) with G{r) oc r~''~^ exp(— r/^). Notice that BCT's work carried 
out for a value of weak anisotropy strength D/J = 3.5 provides a strong confirmation 
for the Itakura's prediction of the coexistence. However, when exploring the nature of 
the low-temperature behavior due to the phases coexistence, questions may arise as to (i) 
whether the memory effects can be "observed" by these simulations and (ii) whether a SG 
transition temperature Tg does exist. In this paper, we only address ourselves to the latter 
question. We find the existence of Tg from a critical slowing-down scaling analysis of the ac 
susceptibility, Xac{T,uj), accompanying some dynamic and ordering features characteristic 
of the low-temperature phases for systems of weak anisotropy. In addition, the behavior of 
the in-phase part of the Xac{T, 00) at a fixed u with various values of D/J provides clues to 
locate the boundaries between anisotropy regions, supporting the Itakura's phase diagram. 

In this work, we perform a MC simulation by using the Metropolis technique {3, Q] , as 
in the BCT's work, for the case D/J = 3.5 in a cubic lattice of L x L x L [L = 10) 
Heisenberg spins. We calculate the temperature dependence of the dc susceptibility in both 
the zero-field-cooled (ZFC), Xzfc(^)) and the field-cooled (FC), xfc(^)) protocols, and of 
the ac susceptibility, Xac(T, tu). At each temperature, an initial A^eq = 100 000-Monte-Carlo- 
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step (MCS) run is carried out to equilibrate the system thermally, followed by another 
Ai'av = 500 000— MCS run to get the time average in each replica. The calculated data are 
then configurationally averaged over a large number of replicas, iVR. The Xzfc(^) and the 
Xfc(^) are calculated in each field protocol by using 



XZFC,FciT) 



1 

H 



(2) 



where H is the strength of the dc field in the range 0.05 < H/J < 0.35, and Mn{t) = 
(S^^/ti • z)/N is the 5-component of the magnetic moment averaged over the t-th spin 
configuration of all N = spins after the equilibrating period (at time < t < A^^av) for the 
R-th replica. To simulate the ac susceptibility, we apply a weak ac field i/ac = HQsm{ujt)z. 
The components of the ac susceptibility, Xac(T,u) = x'{T,uj) — ix"(T, cu), are calculated 
using the formulae in Ref. 9 as 



x\T,co) = — 
-no 



x-"(r,u;) = -— 

-no 



^ _R=1 \ t=l 



COs(co't) 



where Hq/J = 0.05, and the frequency uj is varied in the range lxlO^<u;<lx 
10-1 (MCS-i). 

Figure 1 presents the dc magnetic "measurements" with the ZFC and the PC protocols 
for H/J = 0.08, 0.16, and 0.32. As one can see in Fig. 1(a), the dc susceptibility curves 
resemble a paramagnetic (PM) state in the high-temperature region. A least-squared fit to 
1/Xfc{T) in this temperature region indicates a Curie- Weiss-like transition at Tc/J ~ 1.50, 
as shown in the inset of Fig. 1(a). This behavior causes one to think of the transition 
to a low-temperature ferromagnetic (FM) phase. However, the Curie- Weiss law is ruled 
out in the case of RAM magnetism because a FM long-range order is prohibited in the 
;hree-dimensional RAM model for any arbitrarily small value of anisotropy (i.e., D/J^O) 
171 . Il8[ |. According to Coey, one may refer the ordering phase, at sufficiently low tempera- 
tures, in the RAM to a noncollinear magnetic structure to distinguish it from the collinear 
magnetic structure of the FM phase y. Similar examples for which the Curie- Weis-law 
behavior cannot warrant a low-temperature FM long-range order are ensembles of strongly- 



interacting nanoparticles 
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221]. Also, the concept of domains becomes less obvious 



and unclear in RAM mangetism. For the low D/J region, the domains (ASPEM domains) 



become extremely small for nearest-neighbor interactions. The concept of domains is lost 
altogether in SPEM structures (i.e., in strong D/J region) because the direction of a spin 
changes essentially randomly from site to site, so the domain is reduced to a single site |4|. 
In the weak region, however, one might more or less preserve the concept of a "tiny domain" 
by using a possible hypothesis that at a temperature below Tc, the system consists entirely 
of small magnetic "correlated clusters" or "droplets" which are formed due to FM exchange 
coupling between nearest-neighbor spins. This droplet hypothesis is similar to that in spin 
glasses, where spin-glass domains can be renormalized to give some insight into spin-glass 
dynamics [l2, 16]. As reported in BCT's work 15], the aging effect observed in the current 



case of D/J = 3.5 for the two-time autocorrelation function is a clear evidence strongly 
supporting the droplet hypothesis because the waiting-time dependence and the subaging 



scaling behavior are natural consequences of the droplet theory in spin-glasses [1^, [161, 123 ]. 
As observed in Fig. 1(a), the clusters are probably nucleated somewhere above Tc, where 
the dc susceptibility curves for different applied fields branch (T/J ^ 2.0). This nucleating 
phenomenon may be described by using two competitive processes: the tendency to freeze 
individual spins in random orientations with decreasing temperature (hindering cluster nu- 
cleation), and the alignment of neigboring spins to form correlated clusters by exchange 
coupling. Therefore, H~^{dMpc{T) / dT) may reflect the growth rate of the clusters. As 
shown in Fig. 1(b), the growth rate may reach its maximum, similar to H~^\dMpc{T)/dT\, 
at the corresponding inflection point in the Mpc{T) curve near Tc, which slightly depends 
on H. The existence of SRO clusters near Tc and the history-dependent nonequilibirium 
dynamics (e.g., aging and memory effects) both below and above the peak of Y zfc(T) have 



been strongly evidenced in Yo.7Cao.3Mn03 and Lao.sSro.sCoOa compounds [2J, |25|]. There- 
fore, Tc should be spoken of as a "smeared" transition temperature rather than as the Curie 
temperature of a PM-to-FM transition. 

In the droplet theory, the response of clusters to the field can be considered analogous 
to the thermal activation of NPs in superparamagnets (SPMs). However, in contrast to 
the SPMs the distribution of cluster sizes, -P(t; V"ciuster(i^)), is a dynamic function of their 
time-dependent linear sizes L{t) ~ Kiuster(^)' ^^"^ cluster redistribution at each temperature 
and instant in time leads to a huge number of random energy states for droplet excitations. 
Therefore, some features characteristic of the droplet hypothesis at low temperatures can 
be distinguished from the Xzfc(^) and the xfc(T) data for spin-glass-like systems [23]. 
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As seen in Fig. 1(a), the temperature Tf at which Xzfc(^) and xfc(^) emerge with de- 
creasing temperature almost coincides with the temperature for the Xzfc(^) peak. Another 
interesting observation in Fig. 1(a) is that Xfc(^) displays a clear cusp below Tf and de- 
creases, together with a flat slow decrease of Xzfc(^)) as the temperature is decreased. At 
sufficiently low temperatures, Xfc(^) seems to be nearly constant. The same behavior of 
the dc susceptibility for conventional spin-glasses has recently been understood under the 
random energy model, a phase-space companion of the droplet theory, which was used to 



qualitatively interpret experimental results for a superspin-glass of FesN NPs |23[ |. 

The most interesting phenomenon for the D/J = 3.5 RAM that can be seen in Fig. 
2(a) is the occurence of a double maximum structure in the in-phase part of the ac sus- 
ceptibility, x'iT.oj), while the out-of-phase part, x"iT,uj), displays only one peak. At first 
glance, the shift to low temperature of the high-temperature dispersive peak, T\,{ijj), with 
decreasing frequency seemingly resembles a SPM system of NPs, whose magnetic moments 
have the relaxation times that could be described by an Arrhenius behavior of Neel-Brown 
expression, r = tq exp(£'b/T). However, at temperatures not far below Tc, x'iT.uj) keeps 
on increasing with decreasing temperature to result in a dispersionless peak at To/ J ~ 1.2. 
This dispersionless peak distinguishes the behavior of the ac susceptibility of the RAM from 
that of an ensemble of noninteracting NPs [26] and that of an ensemble of dipolar-interacting 
NPs [9^. As shown in the inset of Fig. 2(b), the so-called Arrhenius plot of J/T\^{uj) vs. 
log(l/u;) indicates a non- Arrhenius behavior with an increasing slope. Notice that the slope 
in a regular Arrhenius plot (i.e., a straightline dependence) provides the value of the activa- 
tion energy that spins need to surmount the energy barrier. Again, under the spirit of the 
droplet hypothesis, the above increase of the activation energy implies a nucleation of small 
clusters below Tq. As a result, the spins in each cluster are so correlated as to hinder free 
thermal activation of individual spins. The growth of clusters with decreasing temperature 
is probably a plausible interpretation for the increasing average energy barrier because a 
strong correlation among spins slows down single-spin thermal activation. Interestingly, the 
conventional critical slowing-down scaling 

can work well with the frequency dependence of T\y{uj). A least-square fit of Eq. (4) to the 
simulated data yields the critical temperature Tg/J ^ 1.20, the critical exponent zv ~ 1.65, 



and the timescale constant r* ^ 10^ (MCS), as shown in Fig. 2(b). To fit Eq. (4) with 
the simulated data, the critical correlation time is obtained by Tc = l/cj, and the reduced 
temperature, e = T/Tg — 1, is produced at each uj when T takes the value of Tb(co'). It is 
shown that the derived value of Tg by the fitting is fairly consistent with its value at the 
dispersionless peak of x'(T,u!). We believe that Tg is the transition temperature to the low- 
temperature CSG phase in which, for weak random anisotropy, there emerges a magnetic 
quasi- long-range order (QLRO) as revealed in Itakura's work 20|], albeit the scaled data of 
x"(T,u!) for e > do not overlap to the universal scaling law, x" (T , u) / XcqiT) = e^GiuTc). 
The same situation has been experimentally reported for an ensemble of interacting Fe-C NPs 
with the concentration of 5 vol% , in which a satisfactory interpretation for the deviation 
from the scaling law was based on the cluster hypothesis. The quite small value of zi/ ~ 1.65 
is attributed to the effect of dynamic finite-size scaling on the value z = zl m such a way 
that the finite-size correlation time should be written as Tc{€,L) = L^^!F{eL^/'^) cx e"^^^ 
rather than as in Eq. (4), which is used for real materials of infinite sizes. For this reason, 
Wansleben and Landau have reported that zlv 1.28 near the critical temperature of a 



cubic Ising ferromagnet [28[, this value is quite comparable to our value. 

To understand whether the dispersionless peak at Tg in x'iT-, ^) is responsible for the onset 
of the transition to the CSG phase in which the magnetic QLRO emerges, we study ensembles 
of nanoparticles with various values of D/ J. As shown in Fig. 3(a), the x'(T, cj) curve at a 
fixed frequency, namely at = 3 x 10^'^ (MCS^^), exhibits the dispersionless peak only for 
anisotropy strength in the range 1.0 < D/ J < 5.0, i.e., in the weak anisotropy region. When 
this region is exceeded, the peak is totally suspended. Strikingly, Itakura, as mentioned 
above, also found a QLRO ground state with a power-law correlation function in this weak 
anisotropy region [2^ . QLRO was theoretically predicted to be a rare case of ordering which 
only exists in weak random-field and random-anisotropy systems of continuous symmetry 
(i.e., systems of non-Abelian symmetry) 29|. The occurence of the dispersionless peak is very 
likely due to the development of the QLRO at temperatures below Tg , which is characterized 
by a smooth rotation of the magnetization over the volume so that the directions of the 
magnetization are ferromagnetically correlated on a quite large or even infinite length, in 
coexistence with a magnetic glassy phase. That may be the reason why Chudnovsky et 
al. fioj] have named this phase correlated spin-glass. Actually, we have recently carried out 
a formal Monte Carlo study of magnetic phase transition for a weak anisotropy strength, 
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D/J = 4.0, of the RAM model. It was shown that a second order transition to the QLRO 



is possible 



30|. It is worth noting that a similar behavior to that shown above for the 



D/J = 3.5 case has been experimental 
anisotropy magnet, by Saito et al 



y observed in a-Ho28Fe72 thin film, a weak random 
3l| . The magnet exhibits not only a critical slowing 
down in phase transition but also a drastical broadening of the distribution of the relaxation 
time around Tg. Moreover, the authors claimed that the dispersionless peak is also observed 
for a sample of a-(Ndi_2:Gda;)28Fe72 films with x = 0, but with increasing x, the peak is not 
observed which suggests a similar physical picture to that accounting for what is observed 
in Fig. 3(a). 

In contrast, low-temperature magnetic behavior is different for systems of anisotropy 
strengths exceeding D/ J ^ 5.0 from that shown above for those of weak anisotropy. In Fig. 
3(a), the dispersionless peak is not observed for anisotropy strengths of D/J > 6.0. This 
disappearance may be understood upon a result in Itakura's work [20^ that the spin-spin 
correlation function for these systems can be well expressed as G{r) oc r^'^'^^ exp(— r/^) 
with rj is an effective exponent and ^ being an effective correlation length. This correlation 
function accounts for a short-range order which is locally attained among neighboring spins 
as a result of scattering spin orientation, and the absence of the dispersionless peak manifests 
a global distruction of the QRLO at low temperatures due to strong random anisotropy 
[19]. In addition, one can see clearly in Fig. 3(a) that when increasing D/J the maximum 
of x'{T,uj) broadens and shifts toward the high-temperature region. The shift, however, 
is reversed backward to the low-temperature region when further increasing D/J beyond 
the value D/J ^ 10.0. Accordingly, it is intuitive to roughly separate the observed short- 
range magnetic behavior of the ac susceptibility into two regions which we preliminarily 
call the intermediate anisotropy region 6 < D/J < 10 and the strong anisotropy region 
D/J > 10, respectively. One question of our concern for these anisotropy regions is that 
whether there exhibits a critical slowing down scaling law as described in Eq. (4) so that 
a spin-glass transition temperature Tg to some magnetic glassy phase still does exist. To 
address this question, we perform our study of the ac susceptibility for another case of 
intermediate anisotropy strength D/J= 10 similar to that of the D/J = 3.5 case. Figure 
3(b) presents the temperature dependence of x'iT,u) for the D/J = 10 case at various 
frequencies. Apparently, the dispersionless peak is not observed for all of the frequencies 
while the high-temperature peak still shifts to low temperature with decreasing frequency 
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similar to what we observed in Fig. 2(a) for the D / J = 3.5 case. To our surprise, the plot of 
J/T\,{uj) vs. log(l/ti;), as shown in the inset of Fig. 3(b), indicates a conventional Arrhenius 
law (i.e., a straigthline dependence), which is quite similar to the behavior of an ensemble of 
noninteracting nonoparticles, and no such a critical slowing down law as given by Eq. (4) is 
fulfilled. The absence of the critical slowing down of a magnetic glassy transition, however, 
contrasts to some experimental results, namely that recently observed by Luo et al. 32^ for 
the Dy-based bulk magnetic glass Dy4oAl24Co2oYiiZr5, a strong random anisotropy magnet. 
The reason is likely due to that the frequencies in our simulation are restricted to vary in 
a so narrow range (one order of magnitude) that only a linear relationship of J/T\^{uj) vs. 
log(l/co') is observed. Further result for a wider range of frequency should be provided to 
have a conclusive reasoning. 

In conclusion, our Monte Carlo simulation gives feacture characteristic of magnetic or- 
dering and collective dynamics for the case D/J = 3.5. With the hypothesis of cluster 
formation, the simulated results show a "smeared" SP-to-SRO transition at Tq » 1.5, fol- 
lowed by an onset of magnetic phase transition at Tg » 1.2 below which both the correlated 
spin-glass and magnetic quasi- long-range order coexist. We demonstrated that the features 
for the D/J = 3.5 case are typically characteristic of those of systems possessing weak 
anisotropy strengths in the range 1.0 < D/J < 5.0. When this range is exceeded, the 
QLRO is predicted to be totally suspended while whether a magnetic glassy transition tem- 
perature Tg does exist is still an open question for further investigations. Our study supports 
a recent theoretical point of view that QLRO is as a rare and novel order that distinguishes 
the weak random-anisotropy glasses from those of strong random anisotopy [29i] . 
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FIG. 1: DC magnetic measurements for both the ZFC and the FC protocols with various apphed 
fields for the D/J = 3.5 case, (a) xzFc(^) and Xfc(^); the best linear fit for T/J> 2.0 (inset) 
indicates a Curie- Weiss-like tranistion at Tq/J ^ 1.5. (b) Myc{T)/H and H''^{dMFc(T)/dT). 
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FIG. 2: Simulated results for the temperature dependence of the ac susceptibility at various fre- 
quencies for D/J = 3.5. (a) x'(T, w) and x"iT,uj) (inset), (b) The critical slowing-down scaling 
log(rc = 1/uj) vs. log(e = 1 — T^/Tg) yields Tg ~ 1.20 and zv ~ 1.65 and the Arrhenius plot of 
J/Tb vs. log(rc = 1/w) (inset) indicates a non- Arrhenius behavior. 
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FIG. 3: Clues for boundaries among anisotropy regions, (a) \'{T,uj) at w = 3 x 10^'^ (MCS^"*^) for 
various D/J values, (b) x'iT,uj) at various uj values. The Arrhenius plot in the inset indicates a 
regular Arrhenius behavior for the D/J = 10.0 case. 
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